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Abstract 

We extend the non-relativistic QCD (NRQCD) prediction for the production of heavy quarko- 
nium with low transverse momentum in hadronic collisions by taking into account effects from all 
order soft gluon resummation. Following the Collins-Soper-Sterman formalism, we resum the most 
singular terms in the partonic subprocesses. The theoretical predictions of J/ip and T productions 
are compared to the experimental data from the fixed target experiments (E866) and the collider 
experiments (RHIC, Tevatron, LHC). The associated non-perturbative Sudakov form factor for the 
gluon distributions is found to be different from the previous assumption of rescaling the quark 
form factor by the ratio of color factors. We further comment on the implication of our results 
on determining the color-octet matrix elements associated with the J ftp and T productions in the 
NRQCD factorization formalism. 
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I. INTRODUCTION 



Heavy quarkonium production has been an important topic in strong interaction physics 
where one can apply perturbative QCD calculations and systematically classify the asso- 
ciated non-perturbative physics in heavy quark systems An effective theory, called 
non-relativistic QCD (NRQCD) fy, has been developed to study the productions and de- 
cays of heavy quarkonia, and especially their productions with large transverse momentum 
in hadronic collisions, up to the next-to-leading order (NLO) in QCD interactions. These 
calculations demonstrated that the so-called color-octet mechanism is important to under- 
stand the J/ip and T productions at the Fermilab Tevatron and the CERN Large Hadron 
Collider (LHC) 

However, how large of the associated color-octet matrix elements remains a question. 
This becomes more important in light of recent studies of the heavy quarkonium polariza- 
tion in the collider experiments, where two separate NLO calculations yield very different 
conclusions The main difference lies on the size of the non-perturbative color-octet ma- 
trix elements. In heavy quarkonium production, the differential cross section is written into 
a factorized form, 

da = J2dHQQ[n}+X)(0 H [n}} . (1) 

n 

Here a(QQ[n] + X) is the cross section at the parton level, which represents the production 
of a pair of heavy quark in a fixed color, spin and orbital angular momentum state n and can 
be calculated perturbatively. The long distant matrix element (LDME) (0 [n]) describes 
the transition of the heavy quark pair in the configuration of QQ[n] into the final state 
heavy quarkonium. The LDMEs are process-independent, which can be extracted from 
experimental data or calculated from non-perturbative method, such as the potential model 
for predicting heavy quarkonium states [l0| • In addition, the LDMEs are organized in terms 
of the velocity v expansion in the NRQCD framework. A fixed order perturbative calculation 
is performed in orders of both the strong coupling constant a s and the power of the velocity 
for the associated LDMEs. For example, for J/ip production in pp collisions, the differential 
cross sections depend on the following three color-octet matrix elements: 

(0 J/ ^[ 3 Sf)]) , [O j/ ^SI\) , (0 J/i '[ 3 P 8 j}} , (2) 

which are at the same order in the velocity expansion. M. Butenschoen and B. A. Kniehl 
used a fixed-order calculation to extract the values of these three matrix elements by fitting 
to several sets of transverse momentum (P±) distributions of J/ip produced in hadron-hadron 
collisions 0,0|. In their fits, the experimental data points in small and mediate P± region are 
all included. They found that the extracted values of these three LDMEs strongly depend 
on the selection of the low P± data sets, by imposing a lower limit on the Px_ of the data 
included in the fits, cf. Tables II, III and IV of Ref. On the other hand, K. T. Chao 
et al demonstrated that J/ip production in high P± region is only sensitive to two linear 
combinations of these three matrix elements [3j. In brief, the above three LDMEs cannot 
be precisely determined from comparing the current experimental data with a fixed-order 
perturbative calculation in the framework of NRQCD. Moreover, the values of these matrix 
elements extracted in 0, @] seem to be too large when compared to those extracted from 
the inclusive cross section of J/ip production in e + e~ annihilation [H, @, @] . 
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In this paper, we will extend the NRQCD prediction for heavy quarkonium production 
to low transverse momentum region with a resummation calculation. This will provide 
additional information on the color-octet matrix elements. Furthermore, comparing the dif- 
ferential cross section spectrum at low P^ and high P± could shed light on the underlying 
mechanisms for heavy quarkonium production, and may also improve our understanding of 
the heavy quarkonium polarization. Some of the recent calculations for the production of 



heavy quarkonia with very high P± could be found in Refs. [ll|, [12|. The low transverse 
momentum heavy quarkonium production has its own interest in applying the perturbative 
QCD factorization formalism. In the low P± region (Aqcd "C P± <C M), there are large 
logarithms of the type a™ ln 2m_1 (M 2 /P^) in high-order perturbative calculations. There- 
fore, to obtain a reliable perturbative prediction, we have to resum these large logarithms. 
In this work we follow the Collins-Soper-Sterman (CSS) formalism [l3j|. We will derive the 
relevant resummation coefficients, by comparing the expansion of the resummation formula 
with the fixed order perturbative calculations within the NRQCD formalism. In particular, 
the most singular contributions from the color-octet and 3 Pj 8 ^ channels will be re- 

summed through the CSS formalism. Earlier work on the soft gluon resummation for heavy 
quarkonium production in low P± region has been performed in Ref. [lij in a model similar 
to the color-evaporation model for heavy quarkonium production. In this paper, we will 
calculate the resummation in the NRQCD framework. 

An important ingredient of the resummation calculation, for predicting low transverse 
momentum distribution of heavy quarkonium produced via gluon-gluon fusion processes, is 
the determination of the needed non-perturbative Sudakov form factor. This form factor was 
previously assumed to be related to that of the quark fusion processes by the ratio of color 



factors Ca/Cf = 9/4 [15|. The latter was determined through a global analysis of Drell-Yan 



(lepton) pair productions in hadron collisions |l6Hl9l|. Our results will show this assumption 



does not work for heavy quarkonium production. Future experimental data on top quark 
pair and Higgs boson productions shall provide additional information on determining the 
non-perturbative Sudakov form factors for gluon-gluon fusion processes. 

The paper is organized as follows. In Sec. II, we derive the low P± behavior of heavy 
quarkonium production in NRQCD by taking the limit of Px_ <§C M . At one-loop order, 
the perturbative corrections are shown to contain large logarithms. We resum these large 
logarithms to all orders in Sec. Ill, following the CSS formalism. The relevant coefficients 
are obtained at the next-to-leading-logarithmic (NLL) level. The numerical evaluations are 
carried out in Sec. IV, where we extract the non-perturbative Sudakov form factor for the 
gluon-gluon fusion processes and the values of the associated color-octet matrix elements in 
the NRQCD framework. We conclude our paper in Sec. V. 



II. LOW P ± BEHAVIOR OF FIXED ORDER CALCULATIONS 

By applying the NRQCD factorization formalism, heavy quarkonium production in 
hadron-hadron collisions arises from the partonic processes, 

a + b -> [QQ] + X , (3) 

where a and b stand for the partons from the incoming nucleons. In high energy collisions, it 
is dominantly produced via the gluon-gluon fusion subprocess. At leading order of v and a s , 
(QQ) pairs are produced in the 1 Sq or 3 Pj configuration via gluon-gluon fusion. They can 
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be produced in either color-singlet or color-octet states, i.e., gg — >■ QQ^S^ 8 } or Q(Q[ 3 -Pj' 8 ]- 
Finite P± arises from the real gluon emission processes. As to be shown later, low P± heavy 
quarkonia are dominantly produced via the 1 5'q' 8 and 3 Pj' 8 channels. This is because in these 
channels the initial state gluon radiation contributes to a singular power behavior 1/P]_ in 
the low P± region, whereas all other channels are power suppressed by P±/M in the limit 
of P± <C M. For example, the 2 — ^ 2 subprocesses for J/ip and T productions include the 
following channels, 

gg ^ [QQ^h + g , (4) 

qg^[QQ J /*> r ] 2 + q, (5) 
qq [QQ J/ *' T h + 9 , (6) 

where [QQ J/ ^ T ] i can be in 1 5'q, 3 Sl' 8 , and 3 P]' 8 configurations, while [QQ J// ^' T ] 2 can be in 
1 S$, 3 Sf, and 3 P}' 8 . The differential cross sections for these channels have been previously 



calculated in Refs. [20|, |21|, |23|, |24j . The low P± behavior of these cross sections can be ob- 
tained by proper expansion in powers of P±/M, where P± and M are transverse momentum 
and mass of the heavy quarkonium, respectively. In this expansion, we only keep the leading 
order contribution, and neglect all higher order terms of P±/M. In the limit P± — > 0, only 
some of these cross sections can reproduce the lowest order terms that are proportional to 
1/P±. This case appears when the respective 2 — >■ 1 partonic subprocess, ab — > [QQ J ^' r ], 
exist. Hence, in high energy hadron-hadron collisions, the production cross section of heavy 
quarkonuim can be approximated as follows: 

a(J/^ ovT) = aigg^'S^ + aigg^'P^ + aigg^'P 8 ) 

+ a(^^ X2 ( 1 P 2 1 ))-Br( X2 ( 1 P 2 1 )^ J/^orT + 7), (7) 

where we have ignored the contribution from qq initial state, qq — > 3 Sf, for its relatively 
small parton density. Furthermore, the production channel gg — > xo{ 3 Pq) 1S n °t included in 
our calculation, for its small decay branch ratio (Br) into J/^> or T. 

In the limit of P± <C M, the differential cross section of the gluon-gluon scattering 
process, gg — > + g, can be expressed as 



d<J 'p x <m = MQ [c] )^r J f{x)dxf{x')dx'^ I 



dyd 2 P ± 



l2\2 



2(1- ft „, . . 

(i-6> + ,(1_6) 



- 1 r x 5(1 - 6) + [2 In - 5 8c ) 5(1 - 6)^(1 - CO 



(8) 



where y and Pj_ are rapidity and transverse momentum of heavy quarkonium, respectively. 
£i = Me v /x\fS, £ 2 = Me~ v /x'yS, and \/S is the center-of-mass energy of the hadron- 
hadron collider, fix) and fix') are parton distribution functions (PDFs). represents 

1 1 8 1 8 

S'q' or 6 Py state, and <5 8c = 1 or for color-octet or singlet channel production. a is 
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proportional to the leading order partonic cross sections, and 



0o 



(q [15 ° 8] ; 



MaW <ofPi)\) 



15 M 7 

5 a 2 s 7T 3 

12 M 5 

4a 2 7T 3 

3 M 7 



(o[ 3 p 2 8 )]) . 



(9) 



It is interesting to compare the above expression to that for the color-singlet scalar particle 
(such as the Higgs boson) production. (See, for example, Ref. |25|.) We find that for the 
color-octet channel, there is an additional term with 5(1 — 6)5(1 — £ 2 )> cf. Eq. OH]), which 
generates sub-leading logarithmic contribution in the low P± region. 

In order to perform the resummation calculation at the next-to-leading logarithmic level 
(NLL), we have to Fourier transform the above expression to the impact parameter space, 
and add the virtual diagram contributions. The Fourier transformation into the impact 
parameter 6-space is defined as W(b) = J d 2 Pj_e~ lP± ' b± dy ^2 P± ■ I n the impact parameter 
space, the logarithmic term will yield a soft divergence in terms of 1/e 2 , where e = (4— D)/2 is 
the dimensional regularization parameter in D-dimension space. This soft divergence will be 
cancelled by the virtual diagrams. After adding up the contributions from real emission and 
virtual diagrams, we are left with collinear divergence in term of 1/e. In our calculation, we 
adopt the modified minimal subtraction (MS) scheme to regularize the remaining collinear 
divergences. Finally, the differential cross sections for describing the production of heavy 



quarkonia in the low P± region 



via gg^ 1 Sl*P J 



1.8 



channels, can be expressed in the impact 



parameter space, within the NRQCD formalism, as 



W {1 \b, M 2 



71 



dxdx' f(x)f(x') 



+(6 -> &)] + S(i - 6)<5(i - & 



-I In* 
2 



M 2 b 2 



6 C A 



b 2 M 2 



1 , 4e" 27B 
e jj, z b z 



(b + -&<&) In 1 



(10) 



where b 



Tpnf)/N c , the expressions for Bq can be found in Ref. [2l||, and f(x 



represents gluon PDF. The gluon splitting function is defined as V, 
x(l-x) + 6(x-l)%. 



<j<j 



x 



(!-*)+ 



+ + 



III. ALL ORDER RESUMMATION 

The large double logarithmic corrections in Eq. fflO|) is a generic feature in low P± dif- 
ferential cross section. In order to have a reliable prediction for low P± heavy quarkonium 
production, we perform all order resummation by following the CSS method. In particular, 
we derive the energy evolution equation for the differential cross section. After solving the 



5 



renormalization group equation, to evolve from the low energy scale of 0(l/b) to the high 
energy scale of O(M), we obtain [13 



W(b, M 2 ) = e-W^A/cWjp U Jj^j , 



(11) 



where the Sudakov form factor is 

dyi 1 



S, 



Sud 



C 2 /b 2 ^ 



111 



CfM 2 



A(C ui i)+B(C x ,C 2 ,ii) 



(12) 



Here C\ and C*2 are two free parameters at the order of unity. The functions A and B can 

oo oo 

be expanded perturbatively in powers of a s , with A = Y ^ and B = (^r) 1 - 

Furthermore, in Eq. (TTTj) . 



8=1 



i=i 



u/^ <n\ch M2 fdxdx' (x x C\ A (x 2 , C\ 2 

xf(x,fi)f(x',n) , (13) 

where x± = Me y / \/S and x 2 = Me~ v / \J~S . After choosing C\ = 2e _TB and C 2 = 1, we 
obtain, for the channels gg — > 1 Sj^\ 3 Pj C ': 



A« = C A , B^ = -(b + -6 c8 )C A 



12 2 



(14) 



, which is the same as that in the gg H 



Together with A^ = C A (§ - g 

process [22], the above result contains all the are needed coefficients for performing the 
resummation calculation at the next-to-leading logarithmic level for the heavy quarkonium 
production in hadron collisions. 

It is interesting to note that the numerical values of the Cgg 1 coefficients for the color-octet 



[cl 

channels are not very different from each other. With the given Bq values from Ref. [21 
we obtain 



CW(^ 8 )=3.16, C«( 3 P 8 ) = 3.76, Cg)( 3 P 2 8 ) = 2.80. 



(15) 



As a consequence, the ratio between the color-octet 3 Pj and 1 5'o contributions, after applying 
the results in Eq.()H]), can be approximately written as 



7.1 ( [ 3 P 8 ]) ^ 7 ( [ 3 P 8 ]) 
Ml (0?S*)]) ~ Mq (OfS*)]) ' 



(16) 



where Mq is the heavy quark mass, and Mq = M/2 in NRQCD factorization formalism. 
Here, we have simplified the two P-wave color-octet contributions by applying the heavy 
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quark symmetry (0[ 3 -P|]) = 5(0[ 3 Pq])- Therefore, at the NLO accuracy, the total contri- 
bution from the above three channels are nearly proportional to the linear combination of 
color octet matrix elements: 

^(Q [8] )~ ^(0[^o 8 ]) + ^(0[ 3 n 8 ])j , (17) 

which is the same as the leading order cross section contributions. The above combination 
will enter into the cross section calculation of the low P± heavy quarkonium production. 
This combination is different from that for calculating the high P± distribution j3|, 0, @]- 
This difference mainly comes from the fact that the gluon-gluon-heavy-quark-pair coupling 
changes with the kinematics. At low P±, we can take heavy quark mass limit (M P±), 
where two gluons are almost on-shell. At large P±, one of the gluon is far off-shell (at order 
of P±) in the 2 — > 2 subprocesses, and the coupling will be different from that in the low 
P± region. These changes depend on the configuration of the heavy quark pair as well. 
Higher order corrections will also change the coupling in the above regions. However, in the 
low P± region, the soft (and collinear) gluon radiation dominates, and the differential cross 
section yields a result proportional to the Born level contribution, cf. Eq.flTT]). This has 
been confirmed by the exact NLO calculation which includes 2—7-3 subprocesses jl, 0]. 



IV. NUMERICAL CALCULATION 

In the numerical calculation, the differential cross section can be written as: 

da = da resum + (da pert - da asym ) , (18) 

where (J IS the resummation of terms which are proportional to 1/P]_ in each order of 
perturbative calculation, and a asym is constructed to cancel the same terms in the fixed- 
order cross section a pert . Therefore, the second term (da pert — da asym ) is power suppressed 
by P±/M in low P± region. The above expression is in principle valid in the whole transverse 
momentum region. In this paper, we focus on the low P± region, where we can safely neglect 
the contribution from the second term. 

In the CSS resummation formalism, the resummation part of total cross section can be 
written as 

■g^W<M = J^'w(b,M, Xl , X2 ) , (i9) 

where W(b, M, xi, x 2 ) has been extensively discussed in the last section. The b integral 
contains contribution from the non-perturbative region where b is so large that a s (l/b) 
cannot be reliably calculated perturbatively. In order to incorporate this contribution, we 
follow the ^-prescription to add a phenomenological non-perturbative form factor [13j, 

W(b) = W{K)W NP (b) , (20) 

where 6* is defined as 6* = bj a/1 + (b/b max ) 2 . Here, the b* does not exceed the b max , which 
is equivalent to making a cutoff on the variable b at b max . 

Generally, the non-perturbative form factor W NP depends on the flavor of the initial 
state partons. There have been many studies in the literature to extract W NP associated 
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with the initial state quarks by comparing theoretical predictions with experimental data 
of Drell-Yan lepton pair production and Z° and W boson productions in pp (pp) collisions. 
(See, for example, Refs. [l6-18||.) On the other hand, for the gluon initiated processes, 



there has been no precise determination of the corresponding W NP (b) in the literature. For 
example, in Ref. [15(, di-photon productions via gluon-gluon fusion processes in hadron colli- 
sions have been investigated, where the associated W NP {b) factor was assumed to scale with 



the color-factor Ca/Cf = 9/4, with respect to that for the quark initiated processes [19 
This is because at the Tevatron energy, the production rate of di-photon events is domi- 
nated by quark initiated subprocesses, so that it is difficult to use that data to extract the 
non-perturbative factor W NP (b) associated with the gluon initiated subprocesses. In this 
work, we will take the BLNY parameterization form, a 3-parameter pure Gaussian form, as 



proposed in Refs. [17H19 



W NP (b) = exp 



-g\ - 5-2 In ( 7j5- ) - gm In (100xiX 2 ) 



(21) 



and extract the non-perturbative form factor associated with the gluon initiated processes 
by performing a global fit to the low P± distributions of heavy quarkonia produced in high 
energy hadron-hadron collisions. More specifically, we will follow the previous studies and 
take the following fixed parameters: Q = 1.6 GeV and b max = 0.5 GeV -1 , in addition to 
the three free parameters g\, gi y and g^. Furthermore, for simplicity, we will assume the 
same non-perturbative function for calculating both the color-octet and color-single heavy 
quarkonium state productions. This is a good approximation because the low P± heavy 
quarkonia are dominantly produced at high energy colliders via color-octet channel. 

As we discussed in the last section, heavy quarkonium (charmonium and bottomonium) 
production rate depends on the associated color-octet matrix elements in the NRQCD fac- 
torization formalism. In particular, we found that the differential cross section in the low P± 
region is proportional to the combination of the two hadronic matrix elements, as shown in 
Eq. (|17p . Hence, we shall treat the very combination of Q} ls o~i and Q[ 3p o 8 l as one free param- 
eter to be determined by the global fits to the J/ip and T production data, separately. We 
note that the hadronic matrix elements are different for J/ip and T productions. In our fits, 
we include the experimental data on the T production from fixed target experiment by the 
E866 Collaboration (26| , and the collider experiments on the T production at the LHC [17 



29] and the Tevatron [30], and J/ifj production at the RHIC [3l| with 1.2 < \y\ < 2.4. The 
other color-sing let matrix elements, (O j ^[ 3 P^)]) and (C T [ 3 P 2 1 )]), are taken from Ref. 
and Ref. [§[ for the J/ip and T production, respectively. In total, our fits contain ten free 
parameters: #1,2,3, two linear combinations of color-octet matrix elements for J/ip and T 
production, respectively, and five normalization parameters Nfu which were introduced to 
account for the normalization uncertainty in each experiment, following the same procedure 
as that in Ref. pi}. 

To determine the free parameters #1,2,3 and the two linear combination of color-octet 
matrix elements in the global \ 2 fits, we minimize the \ 2 contribution from the five experi- 
mental data sets, as shown in Table I, while allowing the normalization of each experimental 
data set to float within one standard deviation of the published experimental error of collider 
luminosity. The result of our fits are shown in Figs. 1, 2 and 3, in which we have plotted the 
shifted theory prediction (multiplied by Nfu for Figs. 1 and 2) with the experimental data 
for the above five experimental data sets. In addition, we also compare the theory predic- 
tion, with the fitted theory parameters, to the J/ip production at the RHIC with \y\ < 0.35 
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1 




12 3 4 

p. (GeV) 



FIG. 1: T production at y/S =7 TeV and \y\ < 2 at the LHC (A), at =1.8 TeV and \y\ < 0.4 
at the Tevatron (B), and at y/~S =38.8 GeV and — 1 < xf < 1 for pp collision at E866 (C). The 
data points are from Refs. SHIS] for E866, CMS, and CDF Collaborations, respectively. 



(cf. Fig. 2(B)) and to the T production at the LHCb (cf. Fig. 3). The fitted parameters 
are summarized in Table I. The uncertainties quoted for #1,2,3 are evaluated with the other 
parameters fixed at their values given by the best fit. The overall agreement between the 
theoretical predictions, based on the NRQCD and soft gluon resummation, and the experi- 
mental data is very good, considering the energy span of the experiments, ranging from fixed 
target to high energy collider experiments. In particular, if we compare the experimental 
data from the Tevatron and the LHC with those from the fixed target experiment, cf. Fig. .1, 
we find that there is a strong energy dependence for the Pi-spectrum, which is reasonably 
described by the energy dependence of the non-perturbative form factor in Eq. f l2"Tj) . Similar 
conclusion holds for the J /if) production between RHIC and LHC, as shown in Fig. 2. 

From Figs. 2 and 3, we find that the agreement between theoretical predictions and 
experimental data for J/ijj production is not as good as that for T production. This may 
be because the mass of J /if) is not very large (as compared to Aqcd) and the theoretical 
uncertainties from the non-perturbative factor of Eq. (12 ip become sizable at relative low mass 
region. We would also like to emphasize that the above resummation results are based on 
the NRQCD factorization formalism for heavy quarkonium production. Hence, the resulting 
non-perturbative factor for the gluon-gluon initiated low heavy quarkonim production 
processes , cf. Eq. (12"TI) . may not be the same as that for the gluon-gluon initiated Higgs 
boson production in hadron-hadron collisions. The former is mainly produced via color- 
octet channel, based on the NRQCD factorization formalism, while the latter is produced 
via color-singlet channel. Although the resummed cross section has properly taken into 
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FIG. 2: J/ip production at \/~S = 
with \y\ < 0.35 (B) and 1.2 < |y| 



7 TeV with 1.6 < \y\ < 2.4 at the LHC (A), at y/S =200 GeV 
< 2.4 (C ) for pp collision at the RHIC. The curve C has been 



multiplied by a factor 1/10 to separate it from the others. The data points are from Refs. 
for CMS and PHENIX Collaborations, respectively. 
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account the different effect of multiple soft gluon radiation for producing a color-octet or 
singlet state in perturbative calculation, a different non-perturbative factor might be needed 
to describe the Higgs boson production. Nevertheless, we shall compare our result to the 
assumption made in Ref. [l8[ for Higgs boson production, 

g x = 0.21 x (9/4) = 0.47, g 2 = 0.68 x (9/4) = 1.53, g x x g 3 = -0.29 , (22) 

which were scaled by the ratio of color factors involved in the gluon-gluon versus quark-quark 
fusion processes. From our fit, we find that the absolute values for the two parameters: g^ 
which controls the scale dependence (InQ term), and the combination g x x g 3 which controls 
the energy dependence (x\X2 term), are both smaller than those in Ref. jlSf j . It will be 
interesting to further investigate this important issue and the relevant phenomenological 
consequences. In order to estimate the uncertainties in our fit, we show in Fig. H] the 
contour plots of the uncertainties of g%, gi, and g\ x g 3 . 

Meanwhile, the differential cross sections in Figs. 1-3 also depend on the color-octet 
matrix elements. As discussed in the last section, in our calculations, we can only determine 
the combinations of the color-octet matrix elements in the form of Eq. f fT7|) . From the fit, 
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Px (GeV) 



FIG. 3: T production at V~S =7 TeV from the LHCb Collaboration. The cures A, B, C, D, and 
E correspond to five rapidity bins: 2 < y < 2.5, 2.5 <y<3, 3<y< 3.5, 3.5 < y < 4, and 
4 < y < 4.5. To separate these cures in one figure, we multiply curves B, C ,D, and E by factors 



1/10, 1/100, 1/1000, and 1/10000, respectively. The data points are from Ref. [29( for the LHCb 
Collaboration. 

we find the following results for J j ip and T productions: 

(0 W [ 1 5 O 8 ]) + ^(0 J/ ^[ 3 P O 8 ]) = 0.0197 ±0.0009 GeV 3 , 



ml 

2 

mi 



<Trl ^o]) + — 2 (0 T [ 3p o]) = 0.0321 ±0.0014 GeV 3 , (23) 



respectively. It's interesting to note that the result of Eq. (T23T) agrees well with that de- 
termined by comparing the total cross section in photoproduction j32j and fixed-target 
hadroproduction [33j of J/ip, which are 0.02 GeV 3 and 0.03 GeV 3 , respectively, as com- 
pared to Eq. fl23|) . Comparing these results with those extracted from the high P± heavy 
quarkonium production 0, [qJ : 



(0 J/ ^[ 1 5 o 8 ])±^(0 J/,/ '[ 3 P o 8 ]) = 0.074 ±0.019 GeV 3 , 

m c 

{0 T [ 1 S^}} + ^4r(0 T [ 3 P o 8 ]> = 0.113 ± 0.020 GeV 3 . (24) 
m b 

We conclude that the color-octet matrix element of [ 3 Pq] is likely to take a negative value. 
Studying the production of heavy quarkonium in high P± region also yields the following 
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FIG. 4: The uncertainty contour for g\ , g\ and g\ x (73. The yellow and blue region correspond to 
Vo CL and 68% CL, respectively 



relation |3|,[9[: 



— n ^6 

(0 W [ 3 S 1 8 ]) + ^P(0 W [ 3 P o 8 ]) = 0.0005 ± 0.00028 GeV 3 



(0 T [ 3 5 1 8 ]) + -^(0 T [ 3 P O 8 ]) = 0.061 ±0.012 GeV 3 . (25) 



Combining equations from Eq. f )23|) to Eq. (125]) . we deduced: 

(0 W [ 1 5 8 ]) = 0.1423 ±0.044 GeV 3 , 
(0 J/,/ fP 8 ])/m;; = -0.0175 ±0.0064 GeV 3 , 

( ^[3 (S -8^ = -0.0093 ±0.0038 GeV 3 , (26) 
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TABLE I: The results of the fits at 68%, 90% and 95% CL. Here, Nfn is the fitted normalization 
factor for each experiment, which is multiplied to the corresponding theoretical prediction to yield 
the comparison presented in Figs. 1 and 2. 



Parameter 


BLNY fit 


gi (CL 68%) 


0.03±0.056 


g 2 (CL 68%) 


0.87±0.065 


g 3 * ff i (CL 68%) 


-0.17±0.011 


9l (CL 90%) 


0.03±0.111 


g 2 (CL 90%) 


0.87±0.080 


53 * 91 (CL 90%) 


-0.17±0.013 


gi (CL 95%) 


0.03±0.113 


<? 2 (CL 95%) 


0.87±0.134 


o 3 * oi (CL 95%) 


-0.17±0.030 


F.866 (T) 

UOUU 1 ± J 


XT — 1 f)« 

J-y fit — i.uu 


(5 points) 


y 2 = 7 9 


Tevatron (T) 


N fit = 0.96 


(5 points) 


X 2 = 3.0 


CMS(T) 


Nftt = 1-05 


(5 points) 


X 2 = 9-2 


RHIC(J/Y>) 


= 0.9 


(6 points) (1.2 < \y\ < 2.4) 


X 2 = 3.4 


CMS(J/V> ) 


N fit = l.H 


(5 points) 


X 2 = 14.6 


x 2 


38 


X 2 /DOF 


1.46 



and 

( T [ 1 1 S O 8 ]) = 0.119 ± 0.021 GeV 3 , 

(0 r [ 3p o]}/ m l = -0.012 ±0.0033 GeV 3 , 

( T [ 3 1 S 1 8 ]) = 0.060 ±0.012 GeV 3 . (27) 

The above result shows that the color-octet matrix elements can be determined to some 
degree by fitting three linear combinations of the LDMEs to the transverse momentum 
distributions of heavy quarkonia (J/tp and T) produced in high energy hadron-hadron col- 
lisions. 

Before closing this section, we note that in the NRQCD factorization formalism, it has 
been assumed that the momentum of the observed quarkonium is the same as the heavy 
quark- ant iquark pair produced at the short distance, where the soft gluon radiation from the 
later stage was neglected. In principle, the (non-perturbative) soft gluon radiation from the 
final state long distance process may affect the P± distribution of heavy quarkonium. It is 
worthwhile to further investigate its effect and implication in heavy quarkonium production 
in hadron-hadron collisions. However, that is beyond the scope of the current paper. 
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V. SUMMARY AND CONCLUSION 



In this paper, we combine the NRQCD and soft gluon resummation formalism to cal- 
culate the P± distribution of J/ip and T production in hadronic collision in the low P± 
region. At high energy colliders, the dominant production mechanism of heavy quarkonium 
is via gluon-gluon fusion, similar to the production of Higgs boson at the CERN Large 
Hadron Collider (LHC). Our analytic calculation shows that the CSS resummation formal- 
ism can be applied to the J /if) and T hadroproduction processes, similar to the Higgs boson 
hadroproduction. Comparing with the color-singlet Higgs hadroproduction, in the CSS re- 
summation formalism, the coefficient has to be modified by adding a constant for the 
color-octet J /if) and T hadroproductions, while the coefficients A^ 1 ' and are the same 
for the color-singlet and color-octet cases. In order to numerically evaluate the Sudakov 
form factor in the large impact parameter (b) region, which is relevant to the low transverse 
momentum (P±) region, we need to introduce a non-perturbative function W NP in the CSS 
resummation formalism. In this work, we use a 3-parameter pure Gaussian form (BLNY 
form) to parameterize the W NP , and fit these free parameters by five experimental data sets 
of J /ip and T hadroproductions, with a total number of 26 data points. We find that we 
need to modify the non-perturbative function W NP previously assumed. The result of our 
analysis is summarized in Table I. Though the non-perturbative function W NP extracted 
from heavy quarkonium data could in principle be applied to studying the distribution of 
Higgs boson produced via gluon-gluon fusion process, it should be emphasized that the main 
contribution of the J/ip and T productions in hadronic collision is from color-octet finial 
state channels, in contrast to the production of the color-single Higgs boson. Hence, they 
might need different non-perturbative factor W NP to describe their low P± distributions. 
Nevertheless, it remains useful to compare this newly determined W NP to that currently 
used for the Higgs boson and di-photon productions in hadronic collisions. 

We have also shown how to extract the values of the color-octet matrix elements from 
studying the transverse momentum distribution of J /if} and T, in both the low and high P± 
regions. The result of this analysis is given in Eqs. ( 126]) and (f2T|) . These matrix elements 
are found to be consistent with the inclusive total cross section in photoproduction [32[ and 
fixed-target hadroproduction [33| of J/if>. Further investigations are needed to clarify the 
underlying mechanisms for heavy quarkonium production in the whole P± range. 

We would like to emphasize that the resummation formula in our calculations are based 
on the extension of the NRQCD factorization in the low transverse momentum region, where 
the most singular contributions are found to follow the CSS resummation expansion at one- 
loop order. It will be interesting to check if this holds at higher orders. Meanwhile, we note 
that the color-single 3 S[^ channel does not contain singular contributions at low transverse 
momentum, because it cannot be produced via a 2 — > 1 process. Hence, we did not consider 
the effect of soft gluon resummation for this production channel in the current analysis. 

VI. ACKNOWLEDGEMENTS 

We thank J. .P. Ma, J. W. Qiu, and in particular, K. T. Chao and his group at Peking 
University for the discussions and comments. We also thank J. C. Peng for discussions on 
the E866 experimental results. This work was partially supported by the U. S. Department 
of Energy via grant DE-AC02-05CH11231, and by the U.S. National Science Foundation 



14 



under Grant No. PHY-0855561. 



[1] N. Brambilla et al. [Quarkonium Working Group Collaboration], |hep-ph/0412 158; N. Bram- 
billa, S. Eidelman, B. K. Heltsley, R. Vogt, G. T. Bodwin, E. Eichten, A. D. Frawley and 
A. B. Meyer et al, Eur. Phys. J. C 71, 1534 (2011). 

[2] G. T. Bodwin, E. Braaten and G. P. Lepage, Phys. Rev. D 51, 1125 (1995) [Erratum-ibid. D 
55, 5853 (1997)]. 

[3] Y. -Q. Ma, K. Wang and K. -T. Chao, Phys. Rev. Lett. 106, 042002 (2011); Phys. Rev. D 
84, 114001 (2011); 

[4] M. Butenschoen and B. A. Kniehl, Phys. Rev. Lett. 106, 022003 (2011); Phys. Rev. Lett. 
107, 232001 (2011). 

[5] M. Butenschoen and B. A. Kniehl, Phys. Rev. Lett. 108, 172002 (2012); K. -T. Chao, Y. - 

Q. Ma, H. -S. Shao, K. Wang and Y. -J. Zhang, Phys. Rev. Lett. 108, 242004 (2012). 
[6] M. Butenschoen and B. A. Kniehl, Phys. Rev. D 84, 051501 (2011) 
[7] M. Butenschoen and B. A. Kniehl, Nucl. Phys. Proc. Suppl. 222-224, 151 (2012) 
[8] E. Braaten, S. Fleming and A. K. Leibovich, Phys. Rev. D 63, 094006 (2001). 
[9] K. Wang, Y. -Q. Ma and K. -T. Chao, Phys. Rev. D 85, 114003 (2012) 

[10] G. T. Bodwin, D. Kang and J. Lee, Phys. Rev. D 74, 014014 (2006) 

[11] Z. -B. Kang, J. -W. Qiu and G. Sterman, Phys. Rev. Lett. 108, 102002 (2012). 

[12] S. Fleming, A. K. Leibovich, T. Mehen and I. Z. Rothstein. far Xiv: 1207.25781 [hep-ph]. 

[13] J. C. Collins, D. E. Soper and G. F. Sterman, Nucl. Phys. B 250, 199 (1985). 

[14] E. L. Berger, J. -w. Qiu and Y. -1. Wang, Phys. Rev. D 71, 034007 (2005). 

[15] see, for example, C. Balazs and C. P. Yuan, Phys. Lett. B 478, 192 (2000); C. Balazs, 
E. L. Berger, P. M. Nadolsky and C. -P. Yuan, Phys. Rev. D 76, 013009 (2007); J. Wang, 
C. S. Li, Z. Li, C. P. Yuan and H. T. Li. larXiv:1205.4311l [hep-ph], 

[16] C. T. H. Davies, B. R. Webber and W. J. Stirling, Nucl. Phys. B 256, 413 (1985). 

[17] G. A. Ladinsky and C. P. Yuan, Phys. Rev. D 50, 4239 (1994). 

[18] F. Landry, R. Brock, G. Ladinsky and C. P. Yuan, Phys. Rev. D 63, 013004 (2001); F. Landry, 
R. Brock, P. M. Nadolsky and C. P. Yuan, Phys. Rev. D 67, 073016 (2003). 

[19] S. Berge, P. M. Nadolsky, F. Olness and C. -P. Yuan, Phys. Rev. D 72, 033015 (2005). 

[20] P. L. Cho and A. K. Leibovich, Phys. Rev. D 53, 150 (1996); Phys. Rev. D 53, 6203 (1996). 

[21] A. Petrelli, M. Cacciari, M. Greco, F. Maltoni and M. L. Mangano, Nucl. Phys. B 514, 245 
(1998). 

[22] Q. -H. Cao, C. -R. Chen, C. Schmidt and C. -P. Yuan. la7Xiv:0909.2305l [hep-ph]. 

[23] M. Klasen, B. A. Kniehl, L. N. Mihaila and M. Steinhauser, Phys. Rev. D 68, 034017 (2003). 

[24] M. M. Meijer, J. Smith and W. L. van Neerven, Phys. Rev. D 77, 034014 (2008). 

[25] X. Ji, J. -P. Ma and F. Yuan, JHEP 0507, 020 (2005). 

[26] L. Y. Zhu et al. [FNAL E866/NuSea Collaboration], Phys. Rev. Lett. 100, 062301 (2008). 

[27] V. Khachatryan et al. [CMS Collaboration], Phys. Rev. D 83, 112004 (2011). 

[28] V. Khachatryan et al. [CMS Collaboration], Eur. Phys. J. C 71, 1575 (2011) 

[29] R. Aaij et al. [The LHCb Collaboration], lar Xiv: 1202.65791 . 

[30] D. Acosta et al. [CDF Collaboration], Phys. Rev. Lett. 88, 161802 (2002). 

[31] A. Adare et al. [PHENIX Collaboration]. larXiv:1105.1966l 

[32] J. Amundson, S. Fleming and I. Maksymyk, Phys. Rev. D 56, 5844 (1997) 



15 



M. Beneke and I. Z. Rothstein, Phys. Rev. D 54, 2005 (1996) [Erratum-ibid. D 54, 7082 
(1996)] 



16 



